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We present the first lattice QCD calculation of the leading-order momentum-independent parity 
violating coupling between pions and nucleons, /i^jviv- The calculation performs measurements on 
dynamical anisotropic clover gauge configurations, with a spatial extent of L ~ 2.5 fm, a spatial 
lattice spacing of Os ~ 0.123 fm, and a pion mass of ~ 389 MeV. While this first calculation 
does not include non-perturbative renormalization of the bare parity-violating operators, a chiral 
extrapolation to the physical pion mass, or contributions from disconnected (quark-loop) diagrams, 
these are expected to result in systematic errors within the quoted statistical error. We find a 
contribution from the 'connected' diagrams of hl'^°^ = (1.099 ±0.5051^; Qg4) X 10 ^, consistent with 
current experimental bounds and previous model-dependent theoretical predictions. 



Quantum chromodynamics (QCD) is the fundamen- 
tal field theory that describes the dynamics and inter- 
actions of quarks and gluons, and the combination of 
QCD and electroweak interactions underlies all of nu- 
clear physics. However, a quantitative understanding of 
nuclear observables directly from QCD has proved elusive 
due to the nonperturbative nature of the theory at low 
energies. Lattice QCD remains the sole avenue for the- 
oretical explorations of observables in the nonperturba- 
tive regime with quantifiable errors. This is particularly 
meaningful for processes which are poorly understood ex- 
perimentally, such as the neutral current parity violating 
(PV) weak interaction between quarks, which is the least 
understood portion of the standard model. In this work 
we report on the first calculation directly from QCD of 
the leading-order momentum-independent parity violat- 
ing coupling between pions and nucleons, ft-^^jv' using 
rif — 2 + 1 lattice QCD calculations on configurations 
with a pion mass of ^ 389 MeV. 

Parity violating interactions have been known since 
the late 1950s[l-3], and their discovery radically changed 
perceptions of the role of fundamental symmetries in par- 
ticle physics. While these interactions can be studied 
in flavor-changing decays, the effects of the PV neutral- 
current in such decays are tiny as the tree-level cou- 
pling between quarks and the Z boson are flavor diago- 
nal and radiative corrections are suppressed by the GIM 
mechanism 0, This leaves PV flavor conserving in- 
teractions as the only laboratories for studying the weak 
neutral current, with the nucleon-nucleon (NN) PV in- 
teraction as the only accessible case. Isolation of the 
hadronic weak neutral current occurs in the A/ — 1 NN 
channel, and this component is thought to be dominated 
by long-range pion exchanged, @. 

At hadronic scales the weak interaction can be consid- 
ered as a pointlike four-quark interaction which gives rise 
to a pion that mediates long-range interactions. Exper- 
iments to uncover this effect are technically demanding 
however, as the ratio of the weak to strong contribu- 
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FIG. 1: (Color online) Model estimates [7l-[l3| (solid line and 
triangles at top) and experimental results (dashed lines with 
labels and la error ellipse in grey, from Refs. [3, fisl . IT^ and 
references therein) for h]^NN versus the dominant isoscalar 
PV coupling combination, along with the results of this work 
(solid vertical line and error band). 

tions to the NN interaction is approximately 10~^. In the 
decades since the discovery of parity violation, a heroic 
series of experiments (see Refs. 0, [l^, [3 ^^'^ references 
therein) have sought to uncover the value of /i^jvj^' 
fined in modern effective field theory language by[6[ 

Z:?.^^ = hl^^{pn+n~nn-p) (1) 

with a proton field p, neutron field n, and pion fields 
7r+/7r^. The most precise of these experiments are plot- 
ted with dashed lines in Fig. [l] with the combined la 
error ellipse shown in grey. The coupling h^j^j^ domi- 
nates the long range parity violating NN potential as it 
is not suppressed by powers of momentum. Although 
lacking precision, experimental results thus far suggest 
that while the isoscalar PV interaction is of natural size, 
the isovector interaction /i^jvAf suppressed. Early re- 
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suits from the most recent experimental collaboration 
to examine nuclear parity violation, the NPDGamma 
collaboration[l3|, have thus far not provided any sig- 



nificant constraint on hi 



NN- 



However, the experiment 



is currently being reinstalled at the Spallation Neutron 
Source at Oak Ridge National Laboratory and should 
soon be able to reach its design precision. 

Because QCD is nonperturbative, how the PV four- 
quark interactions build up into the composite interac- 
tions of the hadrons is not analytically known. Several 
mo del- dependent attempts have been made to calculate 
^ttNN such a way that the nonperturbative effects are 
included. The earliest of these used the quark model 
and symmetry considerations to make the first theoret- 
ical predictions of h^j^j^fl^ (the DDH result). Despite 
tremendous effort, the remaining systematic uncertain- 
ties from the nonperturbative sector of QCD prevented 
Ref. [14] from specifying a result, and instead the out- 
come of the calculation was presented as a 'best guess' 
with an accompanying range of values. Subsequent calcu- 
lations using the quark model[^, , chiral solitons[3, [1], 
and QCD sum rules j7l| - [l3| have obtained greatly varying 
values of h^j^pf, but all have remained within the origi- 
nal DDH range. The DDH range and the results of each 
model calculation are shown at the top of Fig. [T] 

The lattice QCD calculation presented here uses 
anisotropic clover gauge configurations with 2 light quark 
fiavors (the u and d quarks in the isospin limit, m„ = md) 
and one heavier quark fiavor (the s quark) , at a pion mass 
of 389 MeV, spatial lattice spacing of 0.123 fm, and a 
temporal lattice spacing of 0.035 fm[l8|:[i9|. The lattices 
have total dimensions of (2.5 fm)^ x 9 fm. Three-point 
correlation functions of the form 

C'I^bM = {0\OBAt)O%=\t')Ol,{0m (2) 

are constructed, with t the sink timeslice and t' the op- 
erator insertion timeslice. In Eq. [21 the interpolating 
operator 0\ ^ (Obj) is used to create the initial state A 
(destroy final state B) with the quantum numbers of ei- 
ther the proton or the neutron-pion. The proton operator 
is €°-'"^Ua{d^C^5Uc), with color indices a, 6, c. Similarly 
the operator e°'^'^^^Un{d T C'y^U r) creates a neutron-pion 
state (titt) in an S-wave |20l423| . 

Using a three-quark interpolating operator to create 
the nTT state greatly simplifies the contractions, and re- 
moves the need to calculate expensive quark-loop con- 
tributions at the sink which would arise from separate 
n and tt operators. For large Euclidean times t' and 
t — t' the higher energy states induced by these inter- 
polating operators will decay away (as determined from 
analysis of the two point functions of the form CA,B{i) = 
(0|OA,B(t)CA,B(0)|0)), leaving only the proton or nvr S- 
wave state desired. Sandwiched between these operators 
in Eq. [2] is the four-quark operator for the A/ — 1 PV 
interaction. A ratio of the three-point functions to a 
combination of the two-point functions will plateau to 
the constant value of the desired parity violating matrix 
element. 



The four-quark A/ = 1 PV operator can be con- 
structed directly from the standard electroweak interac- 
tion LagrangianQ at the scale of the weak gauge bosons 
by integrating out the Z boson (the contributions from 
the exchange of the bosons are neglected as they are 
suppressed by sin^(6'c) ~ 0.05, where 9c is the Cabibbo 
angle). One can then use continuum one- loop QCD per- 
turbation theory to run the operator coefficients to the 
scale of the hadronic interactions (A^ = 1 GeV) integrat- 
ing out the heavier b- and c-quarks along the wav[6l[24|. 
During the course of this running mixing between opera- 
tors with the same quantum numbers will occur, leaving 
a total of 8 operators at the hadronic scale. There is no 
mixing with lower-dimension operators as the A/ = 1 PV 
operator also conserves CP, precluding quark bilinear op- 
erators from contributing with divergent inverse powers 
of the lattice spacing. The full four-quark A/ — 1 PV 
operator at the hadronic scale can then be expressed as 



-,A/=i _ Gpsm 
'pv — 



1=1 



3V2 



(3) 



where Gp = 1.16637 x 10~^ GeV^^ is the Fermi coupling 
and s\T?{6w) = 0.231 is the weak mixing anglefl^. The 
four-quark operators that contain only light (u and d) 
quarks are 6*^, while the Of contain s-quarks along with 
light quarks. The coefficients Ci and Si of these operators 
and the specific operator forms used for 9i in this work 
can be found in Ref. (26j . 

Performing the quark contractions in the three-point 
correlation function of Eq. [2] using the above operators, 
one arrives at three possible diagrams for the quark prop- 
agators. The first type connects two of the quarks from 
both the source and sink operators to the weak operator, 
with the third quark going directly between the source 
and sink. This type is drawn in Fig. [2{a) and is called 
the 'connected' case. The second, 'quark-loop' type of 
Fig. [2Kb) contains a quark loop at the weak operator in- 
sertion while connecting only one quark each from the 
source and sink to the weak operator. The final type 
contains a weak operator where all four quarks are con- 
tracted with each other, leading to an entirely 'discon- 
nected' contribution. However, in the isospin limit the 
contributions from this type of diagram will sum to zero, 
saving considerable computational expense. Because the 
interpolating operators consist entirely of light quarks, 
the operators Of will have contributions to both the con- 
nected and quark-loop diagrams, while the operators Of 
will contribute only to the quark-loop diagrams as the 
s-quarks will be required to be contained in the quark 
loop itself. 

Typically 3-point correlation functions are computed 
on the lattice using an efficient technique known as se- 
quential inversion, whereby the quark propagators cal- 
culated from the source to the sink are contracted into 
a new 'source' which is inverted to obtain the propaga- 
tor backwards to the operator insertion. However, this 
technique fails for this calculation both in the case of the 
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FIG. 2: The (a) connected and (b) quark-loop diagrams that 
contract the parity-violating operator with the interpolating 
operators for the source and sink. The filled circle and square 
represent the three-quark interpolating operators used at the 
source and the sink respectively, with one positive party and 
the other negative parity. 

connected diagrams (due to the need for two propagators 
between the operator and the sink) and in the case of the 
quark loop diagrams (as the quark-loop would remain to 
be calculated). Instead, this calculation performs two 
separate quark propagator inversions, one at the source 
and one at the weak operator insertion. This method un- 
fortunately restricts the measurements to a single spatial 
site on the operator timeslice (all spatial sites are sam- 
pled over the course of the calculation), but allows for 
maximum flexibility and computational efficiency (as the 
propagators may be used for both the connected and the 
quark-loop diagrams, and for any of the weak operators). 
With this method, the timeslice on which the weak oper- 
ator is placed {t') must be large enough that the excited 
states of the source operator are exponentially small, and 
for this calculation t' = 24. 

As previously mentioned, to extract the desired matrix 
element a ratio of 3-point and 2-point functions must 
be formed such that in the limit of large t' and t — i' 
contamination from excited states dies off and the ground 
state overlap factors are canceled, allowing the ratio to 
plateau to the value of the matrix element. This ratio is 
given by 

c'iit) \c'i{t-t')ci{t)ci{t')) 

where the smearings i and j used in the 2-point functions 
must match that used for the corresponding state in the 
3-point function in order to have the correct cancelation 
of overlap factors. However, as discussed in Ref. [26] the 
differing energy levels of the proton and the nir states 
will cause an insertion of energy by the weak operator to 
occur, modifying Eq. [l]to 

Cpy^ = hlffff (pTT+n - niT^p) + hEDt {pn^n - fnT~p) 

(5) 

with some unknown coefficient He, making the long-time 
behavior of Eq. 0] 

Rrl.^p ^ - {Knn + [Ep - E,,^] ■ He) (6) 

with Ep and En-n the energy levels of the proton and wk 
states, respectively. However, while the energy injection 



term is present in both the forward {p — > nir) and back- 
ward (nTT — >■ p) interactions, it can be eliminated with an 
antisymmetric combination of Eq. [HI leading to a plateau 
region given by 

H^' - \ {R;Un. - Rt^p) ^ hl^N- (7) 

A total of 100,871 measurements of each of the smear- 
ing combinations of H^^ are performed, where i,j can 
be either point- or shell-smearing. These measurements 
are then blocked on each configuration and bootstrapped. 
One can enhance the plateau region for Eq. [7] by taking 
appropriately normalized linear combinations of the dif- 
ferent smearing combinations, using the matrix-pronvpTj 
method on the bootstrapped ensemble to determine the 
optimal linear combination. This is done for both the 
connected and the quark-loop contractions. In the case 
of the quark-loop diagrams the signal-to-noise ratio re- 
mains far too small to recover any reliable result, and 
we do not attempt to extract a signal. It is expected 
that improvements in both contraction algorithms and 
overall calculation runtime will be needed to overcome 
this difficulty and reliably extract the quark-loop con- 
tribution. For the connected contributions the analysis 
returns the data shown in Fig. |31 revealing not only a 
well defined plateau region, but a robust non-zero con- 
tribution to /l^jvAf- 

lOi I I I [ ' I I ' j I J 
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FIG. 3: (Color online) Lattice results for the contribution 
of connected quark diagrams to /i^jvjvi ^ ^ function of Eu- 
clidean lattice time from the operator insertion. The solid line 
is the fully correlated fit value over the plateau region with 
the grey rectangle the statistical plus fit window systematic 
uncertainty. 

In Fig. [31 a fully correlated minimizing fit to a con- 
stant is performed over the plateau region, with addi- 
tional systematic error due to the choice of plateau re- 
gion determined by shifting the ends of the region ±2 
timeslices. The quoted systematic error is one-half of the 
maximum minus the minimum of these shifted fits. The 
fit result and statistical plus systematic error are shown 
in Fig. [3] with the solid line and grey band. The contri- 
bution of the connected diagrams to h^j^j^ is then found 
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to be 

Knn = (1-099 ± 0.505tO;0^«) x 10^^ (8) 

where the first and second uncertainties are statistical 
and systematic, respectively. The fit result is plotted 
in Fig. [T] as the vertical line and error band, and it is 
consistent with both experimental bounds and previous 
model calculations. 

Nonperturbative renormalization of the bare PV oper- 
ators at the lattice scale and subsequent matching to a 
perturbative scheme is not performed for this first cal- 
culation, though results from other four-quark calcula- 
tions (at similar pion mass and lattice spacing) indi- 
cate that this should affect the result by a value sig- 
nificantly below the quoted statistical error [2^. With 
the clover action, lattice spacing errors are expected to 
be 0{alAQ^jj) ~ 2%, also well below statistical uncer- 
tainty. Because sequential propagators are not used, one 
largely eliminates excited state contamination by choos- 
ing an operator insertion time well into the 2-point cor- 
relation function ground state plateaus (though this nec- 
essarily increases the statistical uncertainty by pushing 
the sink operator further into the baryon noise) . Finally, 
one expects from chiral perturbation theory the next-to- 
leading-order contribution to be a pion loop originating 
at the operator insertion, giving an expected finite vol- 
ume error of 0{{m^f^L^)-^) ~ 7% {U = 132 MeV is 
the pion decay constant). While future calculations must 
also address these sources of systematic error, the un- 
certainties in this work remain dominated by statistical 
uncertainty. 



In conclusion, we have performed the first calculation 
of the quantity h\j^j^ directly from the underlying theory 
of QCD. Our calculation was performed on one ensem- 
ble of anisotropic clover configurations with a pion mass 
of ^ 389 MeV. Future calculations will need to be 
performed at pion masses closer to the physical point, 
with sufficient statistical resolution to extract the con- 
tribution of the quark-loop diagrams (expected to be on 
the order of lO^x more measurements), and include non- 
perturbative renormalization of the included operators. 
While significant technical challenges remain in the cal- 
culation of the full matrix element, this first of its kind 
result clearly shows that lattice QCD can make a signifi- 
cant contribution to the theoretical, model-independent, 
understanding of quantities that are difficult to access 
experimentally. Our initial result shows good agreement 
with current experimental bounds and paves the way to- 
ward a complete extraction of h\j^j^ at a precision con- 
sistent with, or better than, the anticipated results of the 
upcoming NPDGamma experiment at Oak Ridge. 
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